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Abstract 

We provide a comprehensive test of factorization in the heavy-heavy B decays motivated by the 
recent experimental data from BELLE and BABAR cohaborations. The penguin effects are not 
neghgible in the B decays with two pseudoscalar mesons. The direct CP asymmetries are found 
to be a few percent. We give estimates on the weak annihilation contributions by analogy to the 
observed annihilation-dominated processes. The Nc insensitivity of branching ratios indicates that 
the soft final state interactions are not dominant. We also study the polarizations in B D*D*^^^^ 
decays. The power law shows that the transverse perpendicular polarization fraction is small. 
The effects of the heavy quark symmetry breaking caused by the perturbative QCD and power 
corrections on the transverse polarization are also investigated. 
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I. INTRODUCTION 



The study of B meson weak decays is of high interest in heavy flavor physics and CP 
violation. In particular, much attention has been paid to the two-body charmless hadronic 
B decays, but there are relatively less discussions on the decays with charmful mesons, such 
as the modes with two charmed-meson in which the final states are both heavy. However, 
the two charmed-meson decays can provide some valuable and unique information which 
is different from the light meson productions. For example, CP asymmetries in the decays 
oi B ^ /}(*)+/}{*)- play important roles in testing the consistency of the standard model 
(SM) as well as exploring new physics fli]. Moreover, these decays are ideal modes to check 
the factorization hypothesis as the phenomenon of color transparency for the light energetic 
hadron is not applicable. Since the decay branching ratios (BRs), CP asymmetries (CPAs) 
and polarizations of -B D^*'> D^*\D^*'> oi*^) have been partially observed in experiments 



it is timely to examine these heavy-heavy B decays in more detail. 
At the quark level, one concludes that the two charmed-meson decays are dominated by 
tree contributions since the corresponding inclusive modes are b — * qcc with q = s and d. It 
is known that the factorization has been tested to be successful in the usual color-allowed 
processes. However, the mechanism of factorization in heavy-heavy decays is not the same 
as the case of the light hadron productions. The color transparency argument 5] for light 
energetic hadrons is no longer valid to the modes with heavy-heavy final states. The reason 
can be given as follows. Due to the intrinsic soft dynamics in the charmed-meson, non- 
vanishing soft gluon contributions are involved in the strong interactions between an emission 
heavy meson with the remained BD^*^ system. Since the corresponding divergences may 
not be absorbed in the definition of the hadronic form factor or hadron wave function, the 
decoupling of soft divergences is broken. This means that the mechanism of factorization has 
to be beyond the perturbative frameworks, such as QCD factorization jol] and soft-coUinear 
effective theory The large Nc limit is another mechanism to justify factorization j^, 
corresponding to the effective color number A^^c = oo in the naive factorization approach 
Q]. The understanding of factorization in heavy-heavy decays requires some quantitative 
knowledge of non-perturbative physics which is not under control in theory. In this paper, 
we will assume the factorization hypothesis and apply the generalized factorization approach 
(CPA) |lfll. to calculate the hadronic matrix elements. 
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It is known that annihilation contributions and nonfactorizable effects with final state 
interactions (FSIs) play important role during the light meson production in B meson decays. 
For instance, to get large strong phases for CP asymmetries (CPAs) in B Ktt and B vrvr 
decays, these effects are .uc.uded inevitably Q Q. Moreover, they are also crucial to 
explain the anomaly of the polarizations in B ^ (j)K* decays, measured by BABAR [1J| 
and BELLE j3] recently. By the naive analysis in flavor diagrams, one can easily see that 
the decay modes of 5 ^ /)(*)0/)(*)o ^nd D^f^^D^*^^ are annihilation-dominated processes. 
Therefore, measurements of these decays will clearly tell us the importance of annihilation 
contributions in the production of two charmed-meson modes. 

For the color-allowed decays, since the short-distant (SD) nonfactorizable parts are as- 
sociated with the Wilson coefficient (WC) of Ci/Nc, where Ci is induced by the gluon-loop 
and it is much smaller than C2 ~ 1, we can see that the effects arising from the SD non- 
spectator contributions should be small Nevertheless, long-distant (LD) nonfactorizable 
contributions governed by rescattering effects or FSIs may not be negligible. Inspired by the 
anomaly of the large transverse perpendicular polarization, denoted by R±, in B —>■ (pK* 
decays, if there exist significant LD effects, we believe that large values of R± may appear 
in i? ^ D*D* and B D*D* too. As we will discuss, the power-law in the two-vector 
charmed- meson decay leads to a small R±. The recent measurement of the polarization 
fraction by the BELLE collaboration gives R± = 0.19 ± 0.08 ± 0.01 [s^ in which the central 
value is about a factor of three comparing with the model-independent prediction within 
the factorization approach and heavy quark symmetry. Clearly, to get the implication from 
the data, we need a detailed analysis in two charmful final states of B decays. 

To estimate the relevant hadronic effects for two-body decays in the B system, we use 
the GFA, in which the leading effects are factorized parts, while the nonfactorized effects 
are lumped and characterized by the effective number of colors, denoted by A^^*^. Note that 
the scale and scheme dependences in effective WCs Cf^ are insensitive. In the framework 
of the GFA, since the formulas for decay amplitudes are associated with the transition form 
factors, we consider them based on heavy quark effective theory (HQET) We will also 
study their Og Q] and power corrections which break heavy quark symmetry (HQS) jioj ]. 
In our analysis, we will try to find out the relationship between the HQS and its breaking 
effects for R±. In addition, we will reexamine the influence of penguin effects, neglected 
in the literature bol. We will show that sizable CPAs in S° D+D- and B- D^D' 
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decays may rely on FSIs. 

This paper is organized as follows. In Sec. II, we give the effective Hamiltonian for the 
heavy-heavy B decays. The definition of heavy-to-heavy form factors are also introduced. 
In Sec. Ill, we show the general formulas for B to two charmful states in the framework of 
the GFA. The effects of the heavy quark symmetry breaking on the transverse perpendicular 
polarization are investigated. In Sec. IV, we provide the numerical predictions on the BRs, 
direct CPAs and the polarization fractions. Conclusions are given in Sec. V. We collect all 
factorized amplitudes for B PP, PV{VP) and VV decays in Appendixes. 



II. EFFECTIVE INTERACTIONS AND PARAMETRIZATION OF FORM FAC- 
TORS 



The relevant effective Hamiltonian for the B meson decaying to two charmful meson 



states is given by 

H,s{^B = l) = ^{K,\CJC'i(/^)Or + C^2(/i)0^] + K.K;[C^i(/^)0? + C'2(/x)0^] 

10 ^ 

-Vi,v,;5^c,(/i)0fcU//.c. , (1) 

k=3 J 

where q = s and d, Vij denote the Cabibbo-Kobayashi-Masikawa (CKM) matrix elements, 
Ci{fi) are the Wilson coefficients (WCs) and Oj are the four-fermion operators, given by 

Oi = {qiUj)v-A{Ujhi)v-A, O2 = iqiUi)v-AiUjbj)v-A, 

01 = {qiCj)v-A{cjbi)v-A, O2 = {qiCi)v-A{cjbj)v-A, 



03(5) = {qibi)v-A ^{q'jq'jhTA, 04(6) = {qibj)v~A ^{qjq'i}vTA, 

q' q' 

3 3 
^7(9) = i^{qibi)v-A^eq,{q'^q'j)v±A, 08(io) = -{qibj)v-A^eq>{qjq'^v±A- (2) 

q> q> 

with i and j being the color indices, O^^q {Ot^iq) the QCD (electroweak) penguin operators 
and {qiq2)v±A = qilfii^ ± 75)92- In order to cancel the renormalization scale and scheme 
dependence in the WCs of Cj(/i), the effective WCs are introduced by 

C(/i)(0(/i)) = C^fif(0)tree- (3) 
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Since the matrix element (O)tree is given at tree level, the effective WCs are renormalization 
scale and scheme independent. To be more useful, we can define the effective WCs as 



eflf _ ^eff I '^1 ^eff _ ^eff , '^2 

eff(q) _ ^eff , '^^3 , 3 / ^eff , 



2 

eff(q) _ ^eff , ''^4,3 O I r^S , ''^lO.g 

eff(q) _ ^eff , ^6,5 . ^ I ^ , "-^S,? 



where 
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with Xi being the non-factorizable terms. In the GFA, l/{Nf^)i are assumed to be universal 
and real in the absence of FSIs. In the naive factorization, all effective WCs Cf^ are reduced 
to the corresponding WCs of Cj in the effective Hamiltonian and the non-factoziable terms 
are neglected, i.e., Xi = 0- 

Under the factorization hypothesis, the tree level hadronic matrix element (O)tree is fac- 
torized into a product of two matrix elements of single currents, which are represented by 
the decay constant and form factors. The B ^ H [H — D, D*) transition form factors are 
crucial ingredients in the GFA for the heavy-heavy decays. To obtain the transition elements 
of S — > if with various weak vertices, we first parameterize them in terms of the relevant 
form factors under the conventional forms as follows: 

{D\V,\B) = F,{q^)[p, - ^q,] + ^Fo{q')q„ (6) 

TUB 'T TUjj* 



{D*{e)\A,\B) 



2mD*Ao{q^)—^q^ + {uib + m£,.)^i(g^) (^e* ^^^j 

ruB + tud* ^ q^ 



(7) 



where V^, = q-ff.b, = q^/)^h, mB,D,D' are the meson masses, denotes the polarization 
vector of the D* meson, P = ps +Pf)(*), q = Ps — PdM and P ■ q = m\ — "m^j,). According 
to the HQET, it will be more convenient to define the form factors in terms of the velocity 
of the heavy quark rather than the momentum. The definition of these form factors can be 
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found in Ref. and the relation with the conventional ones are given by 

2, 1713 + rUD 



A,{q') 



2y/mBmD 
rriB + tud 

rriB + niD- 



niB-mD , . 
^+(^) - - — — ^Aw) 



Coiq') 



rriB + tud 

TUB + TTLd f W — 1 



, Mq 



TUB — mo \w + 1 

2\ _ + 



2y/mBmD* 
niB + rriD* 
2y/rnBmo* 

■rriB + mo* 

2y/mBmD' \ niB + mo- 
rriB - mo 



2^fm^mY)* 



niB 



rriB 



2m 



D* 



Ao{q') = A,{q' 



q 



mB 



mB 
mp* 
mB 



(8) 



AmBmo* V 

withcij = {m\+m\j — q^) / {2mBmH) andCH(g^) = l—q^/{mB+mH)'^- It is known that under 
the HQS, ^+ = = ^Ai = ^A-i = ii.'U)) and ^_ = = 0. In our numerical estimations, we 
will base on the results of the HQS and include and 1/mB power corrections as well. 



III. GENERALIZED FACTORIZATION FORMULAS AND POLARIZATION 
FRACTIONS OF VV MODES 

By the effective interactions and the form factors defined in the previous chapter, the 
decay amplitude could be described by the product of the effective WCs and the hadronic 
matrix elements in the framework of the GFA. For the hadronic matrix elements in i? — > PP 
decays, we will follow the notation of Ref. given by 

j^iBP,,P2) ^ (P2|feg3)v'^A|0)(Pi|(gi6)v'-A|5)=^/p,(m|-m|jFo^^^(m?,J, 

Xr^'^^) ^ {P,\iq,q,)s^p\0){P^\iq^b)s-p\B) = -i^^ f p ^^^ ~ '^^''^ F^^m^ (9) 

m2 + ms mfo — m\ 

where {qib)s-p = gi(l -75)6, {q2qz)s+p = ^2(1 + 75)^3 and mf,,i,2,3 correspond to the masses 
of quarks. The vertex {S — P)®{S+P) is from the Fierz transformation of {V — A)®{y + A). 
To get the decay constant and form factors for scalar vertices, we have utilized equation of 
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motion for on-shell quarks. Moreover, we use 

xiBPy) ^ (^v\{q2q3)v^A\0){P\{qib)v-A\B) = 2fymvFf''{ml.){e* -p,), 
X(^^-^) ^ (P|feg3)y-A|0)(r|(gi6)y-A|5) = 2/pmv<^(m|)(e*-pJ, 
Xf^'^) - {nq2q3)s+p\0){V\ms-p\B) = ^^fp^^A^imlW .p,), (10) 

1712 + ^3 n^b + ^1 

and 

^(ijy.y.) ^ {V2\{q2q3)v-A\0)mq,b)y_A\B) 



2 



+^e^,ya|3£ 2^1 P2P1 



2 



(11) 



for i? PV{VP) and i? ^ respectively. We note that the sign difference of 

^<^Pi,P2) ^(j3Pi,P2) Q ^jjj make the penguin effects become non-neghgible. 

On the other hand, the same sign of x'^^'^^ and x!^^'^^ in Eq. (jiup leads to the 
penguin effects negligible. Since the time-like form factors for annihilation contributions 
are uncertain, we take F/^f^^^^^'^ = {MiM2\{q2qz)vTA\Q) {Q\{qih)v-A\B) and Fg^^'^^^^) ^ 
{MiM2\{q2qz) s+p\^) (0|(^i&)s-p|-B) to represent them, where M can be pseudoscalar or vec- 
tor bosons. Note that due to the identity of ^^(pi) -pi = 0, we have used {V\{q2q3)s+p\^) = 0. 
With these notations and associated effective WCs, one can display the decay amplitude for 
the specific decay mode. We summary the relevant decay amplitudes in Appendixes. Once 
we get the decay amplitude, denoted by A{B M1M2), the corresponding decay rate of 
the two-body mode could be obtained by 

TiB ^ M1M2) = -^^\AiB ^ MiM2)p. (12) 

with p being the spatial momentum of Mi^2- Consequently, the direct CPA is defined by 

^ ^ r(^ -> M1M2) - T{B ^ M1M2) 

r{B M1M2) + T{B M1M2) ■ ^ ^ 

Besides the BRs and CPAs, we can also study the polarizations of the vector mesons in 
B —>■ VV decays. To discuss the polarizations, one can write the general decay amplitudes 
in the helicity basis to be 



(14) 
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In this basis, the amphtudes with various hehcities can be given as 



Hq = —ax — b{x'^ — 1), H± = a ± V — 1 c. 

where x = {m\ — rriy^ — rriy^) / {2'my^mv2) ■ addition, we can define the polarization 
amphtudes as 

= -^{H+ - H^) = V2Vx^ - 1 c. (15) 
v2 

Accordingly, the decay rate expressed by helicity amplitudes for the VV mode can be written 

as 

and the polarization fractions can be defined as 

= I /I 12 , I /l2 , I /I 12- (16) 



where i = and || (_L), representing the longitudinal and transverse parallel (perpendicular) 
components, respectively, with the relation of Rq + R\\ + R± = 1. Under CP parities, i?o,|| 
are CP-even while R± is CP-odd. 

From the hadronic matrix element in Eq. (fTTj) . the amplitudes a, b and c in the framework 
of the GFA are expressed by 



a = -CeffirriB + mvi)mv2fv2Af^^(m 



2mv, ml 



where Ce// represents the involved WCs and CKM matrix elements. With the form factors 
in Eq. (jH)) and the heavy quark limit, we get that the ratios = b/a and = c/a are 
related. Explicitly, we have 

H = r,= . ^0.16, 18 

niB (1 + w) 

which are small. From Eqs. (fT3j) and (fT^ . we find that the polarization fractions behave as 



m 



m 



B 



i?o~i?||, , (19) 



which indicate that the power law in the heavy-heavy decays is different from the hght-hght 
ones, which are expected to be i?o ~ 1, R\\ ~ R± ~ (^(^y/^^s)- Moreover, R± is directly 
related to c and can be written as 



Ra 

with To = 1 + (x^ - l)|rcp + \x + (x^ 
we find an interesting relation 

Ra 



1 



(20) 



2 

l)rh\ /2. By comparing to the result in the HQS, 



R 



HQS 



Cd* 



(21) 



where R^^^ = 0.055 denotes the transverse perpendicular fraction under the HQS. As a 
good approximation, the form factor y42-dependent of R± is decoupled. By the relationship, 
we can clearly understand the influence of the HQS breaking effects. 



IV. NUMERICAL ANALYSIS AND DISCUSSIONS 
A. Estimations on the annihilation contributions 

Since the annihilation contributions relate to time-like form factors and there are no direct 
experimental measurements, we shall neglect them in our calculations. However, to make 
sure that the neglected parts are small, we can connect the processes of i? ^ D^D^ and B — >■ 
D^D^ to the decays DjK~^ and B^ J/ipD^, which are directly associated with 

annihilation topologies, with the experimental data of Br{B^ D^K^) = (3.8±1.3) x 10^^ 
P] and Br{B^ — J/ipD^) < 1.3 x 10~^ js^, respectively. By the flavor- diagram analysis, 
shown in Fig.[TJ except there appears a CKM suppressed factor Vcd ~ A (see FigH^) for D^D^ 



q' q 



(a) 



(b) 



FIG. 1: Flavor diagrams for (a) Bd D+K-{J/^D^) decays while q' = s{c) and (b) B 
D^D^{DfDj) decay while q = u{s). The gluon attached denotes the nonfactorized effects. 
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and D^Dg modes, the four modes have the same decay topologies. Hence, by assuming that 
they have similar hadronic effects, the BRs of B*^ —>■ D^D^[DfD~) could be estimated to 
be less than 0(10-^). 

To give a detailed analysis, we can include the character of each mode, governed by 
the meson distribution amplitudes. For simplicity, we will concentrate on the leading twist 
effects and take the meson wave functions to be $d oc fDx{l — x){l + 0.8(1 — 2x)) j^l, 



$n cx foM^ -x){l + 0.3(1 - 2x)) j2j, J/^ oc fj/^x{l - x){x{l - x)/(l- 2.8x{l - x)))^-^ 
y and oc //^x(l - x)(l + 0.51(1 - 2x) + 0.3[5(1 - 2x)2 - 1]) [2^, where x is the 
momentum fraction of the parton inside the meson and fD,Ds,J/^,K are the decay constants 
of D, Dg, J/\l' and K mesons, respectively. From these wave functions, we know that 
the maximum contributions are from xq ~ (0.35,0.43,0.5,0.5) for {D, Ds, J/"^ , K). With 
the information, we can estimate the decay amplitudes in order of magnitude for Bd — »• 
D+K-{J/^D^) and Bd D'^D^D+D;) as shown in Figs. CJi andCb with q' = s(c) and 
q = u{s), respectively. Note that there exists a chiral suppression in the factorized parts 
in annihilation contributions. However, we just consider the nonfactorized effects in the 
estimations. Therefore, by Fig. Q with the gluon exchange, the decay amplitude is related 
to the propagators of the gluon and the light quark, described by 1/(A;2 + ^3)^/(^2 + ^3)^, 
where ^2(3) denote the momenta carried by the spectators. For simplicity, we have neglected 
the momentum carried by the light quark of the B meson. By the momentum fraction, the 
decay amplitude could satisfy that A oc l/(x2X3)^. Hence, the relative size of the decay 
amplitudes could be given approximately as 

AiDjK-) : AiJ/^D') : A(D°D°) : AiD+D;' - ^^^^^ ■ ■ ■ 



(X2X3)2 ■ (X2X3)2 ' (X2X3)2 ' (X2X3)2 ' 

With the information of maximum contributions, characterized by xq for each mode, we get 



A{DjK-) : A(J/^D°) : A(D°D°) : A{DjD; 

(22) 



^ /d/j^/O^Y Xpn / 0.5-0.43 Y Xfn, / 0.5 
■ JdJk \0.Q5j ■ JdJk [ 0.352 ) ■ (^0.43 



With the kinetic effects, the ratios of BRs are roughly to be BR{D+K~) : BR{J/^!D^) : 
BR{D^D^) : BR{DfD;) ~ 1 : 0.25 : 0.39 : 0.16. That is, the BRs of 5 ^ D^D°{D+D;) 
could be as large as O{10~^), which implies that annihilation effects could be neglected in 
the discussions on the BR of the production for color-allowed two charmful mesons. We note 
that our estimations are just based on SD effects and at the level of order of magnitude. 
Clearly, direct experimental measurements are needed to confirm our results. 
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B. as, power corrections and the parametrization of Isgur-Wise function 



In the HQS limit, the form factors could be related to a single Isgur-Wise function ^{oo) 
by i+ = iv = = iA-i = iiw) and = = 0- We now include the perturbative 
QCD corrections induced by hard gluon vertex corrections of 6 — c transitions and power 
corrections in orders of l/mq with Q = b and c. Consequently, the form factors can be 
written as 

Uw) = {a, + /3^iw) + %iw)} ^(w). (23) 

where ^{w) is the Isgur-Wise function, a+ = ay = aAi = = 1, a- = = and j3i{uj) 
and 7j(c(j) stand for effects of ag and power corrections, respectively. Explicitly, for the two- 
body decays in our study, ~ 1.3 and the values of the other parameters are summarized 
as follows HQ: 

P+ = -0.043, p_ = 0.069, Pv = 0.072, 

(3^^ = -0.067, (3a, = -0.114, (3a, = -0.028, 

(24) 

7+ = 0.015, 7„ = -0.122, 7y = 0.224, 

7a, = 0.027, 7^2 = -0.093, 7A3 = 0.014. 

Clearly, the range of their effects is from few percent to 20% level. In particular, the power 
corrections to the form factor (or say V^q"^)) are the largest, about 20%. The resultant 
is also consistent with other QCD approaches, such as the constitute quark model (CQM) 
and the light-front (LF) QCD 

After taking care of the corrections, the remaining unknown is the Isgur-Wise function. To 
determine it, we adopt a linear parametrization to be C,{w) = 1 — p^(w — 1) for the transition 
B ^ H, where p'jj is called the slope parameter. We use the BRs of semileptonic B 
D^*H~i'i decays to determine p|^. We note that the values of pj^ and p|,, are not the same 
as those in D and D* decays. In our approach, the difference is from higher orders and power 
corrections. Hence, with = 0.0416 and BR{B D^*H-ue) = 2.15±0.22(6.5±0.5)% Q, 
we obtain pj^ = 0.90 ± 0.06 and pj^, = 1.09 ± 0.05. Since the errors of p^^ are small, we will 
only use the central values in our numerical results. 
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C. Results for BRs and polarization fractions 



To get the numerical estimations, the input values for the relevant parameters are taken 



to be as follows 



J 29, 3: 



/b = 0.20, fD*= 0.23, /d. =0.24, /^j = 0.275 GeV; 
Vcd = -A, = 1 - AV2, = AX\ 

Vt, = X\V,b\Rte-"^\ Vts = -AX\ Vtb = l, 
A = 0.83, A = 0.224, 0i = 23.4°, Rt = 0.91; 

ruu = 0.005, rud = 0.01, = 0.15, rUc = 1.5, rub = 4.5 GeV. (25) 

Note that the numerical results are insensitive to light quark masses. As to the WCs, we 
adopt the formulas up to one-loop corrections presented in Ref. jll| and set /i = 2.5 GeV. 
As mentioned early, since the nonf actor ized contributions are grouped into Nf, the color 
number in Eq. will be regarded as a variable. To display their effects, we take the values 
of Nf = 2, 3, 5 and oo. 

By following the factorized formulas shown in Appendixes, we present the BRs with 
various A*"^^ in Tables HI ITTl and UTTl for PP, PV{VP) and VV modes, respectively. In order 
to accord with the experimental data, our predictions of the BRs are given as the CP- 
averaged values. For comparisons, we also calculate the results in terms of the form factors 
given by the CQM and LF, which are displayed in Table IIVI Since the CPAs are quite 
similar in different models, in Table HVI we just show the results in our approach. As to the 
polarization fractions, we present them in Table Therein, to understand the influence of 
the HQS breaking effects, we separate the results to be HQS and HQSi(ii), representing the 
HQS results and those with as (as+power) corrections, respectively. 

TABLE I: BRs (in unit of 10"^) for B ^ PP decays with = 0.90. 





mode 


Nf = 2 


Nf = 3 


Nf = 5 


Nf = oo 


Exp. 




- D+D- 


7.26 


8.25 


9.06 


10.46 


8±3\2] 


B- 


- D'D- 


7.85 


8.94 


9.82 


11.34 


13±4 [2j 


5° 


D+D- 


0.28 


0.31 


0.34 


0.40 


0.321 ±0.057 ±0.048 \3] 


B 




0.30 


0.33 


0.37 


0.43 


0.562 ± 0.082 ± 0.065 [3] 
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TABLE ir 


BRs (in 


1 unit of 10^'^') 


for B PV{VP) decays with /j^j,. 


= 90fl 091 




mode 




2 Nf = 3 


iV^ = 5 


N^ = oo 




Exp. 




D+D*- 


9.52 


10.80 


11.84 


13.62 




10 ± 5 


5° 


D*+D- 


6.78 


7.67 


8.41 


9.66 




10.7 ±2.9 \2] 


B 


D^D*- 


10.35 


11.73 


12.87 


14.79 




9±4 [2] 


B- 


D*^D- 


7.37 


8.34 


9.14 


10.49 




12 ± 5 [2] 


B^ 


D*+D- 


0.25 


0.29 


0.32 


0.36 






B^ 


D+D*^ 


0.37 


0.42 


0.46 


0.53 






5° 


D*+D- 
+ D+D*- 


0.62 


0.71 


0.78 


0.89 




0.93 ±0.15 [2j 


B 


D^D*- 


0.40 


0.45 


0.50 


0.57 


0.459 ± 0.072 ± 0.056 [3] 


B 


D*°D- 


0.28 


0.31 


0.34 


0.39 







TABLE IIL BRs (in unit of lO^^) for B ^ VV decays with pf^. = 1.09. 





mode 


Nf = 2 


Nf = 3 


Nf = 5 


Nf = oo 


Exp. 




^ D*+Dl- 


22.52 


25.51 


27.98 


32.19 


19 ±5 a; 














18.8 ±0.9 ±1.7 f4] 


B- 


D*^D*- 


24.44 


27.69 


30.37 


34.93 


27 ± 10 [2] 


B^ 


D*+D*- 


0.87 


0.91 


0.99 


1.14 


0.81 ±0.08 ±0.11 [3]; 














0.87 ±0.18 [2j 


B 


D*°D*- 


0.81 


0.98 


1.08 


1.24 





We now present our discussions on the results as follows: 

(1) The non-factorizable contributions are not dominant for color-allowed two charmed- 
meson decays. According to the classification in Refs. , the decay modes displaced 
in Tables IH ITTl and IlIII belong to class I, which are dominated by the external ly-emission. 
The leading decay amplitudes are proportional to the effective coefficient ai, which is stable 
against the variation of A^^^^. Thus, the predicted branching ratios are insensitive to Nf . 
This means that annihilation contributions and FSIs, neglected in the GFA, are sub-leading 
contributions. On the other hand, by varying Nf from 3 to 2, or 3 to oo, the branching ratios 
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TABLE IV: BRs (in unit of 10"^) with Nf = 3 and the form factors calculated in the CQM, LF 
and HQSC, while the CPA is only shown in our approach (HQSC). 





Mode 


CQM 


LF 


HQSC 


Acp{%) 




^ D+D- 


9.70 


10.33 


8.25 


-0.2 


B- 


^ D'^D- 


10.58 


11.26 


8.94 


-0.2 




D+D- 


0.38 


0.40 


0.31 


2.5 


B- 


-^D'^D- 


0.42 


0.44 


0.33 


2.5 


50 


D+D*^ 


12.49 


11.42 


10.80 


—0.1 


B^ 


D*+D- 


9.19 


8.50 


7.67 


0.0 


±J 






1 9 47 


1 1 73 


— n 1 


5- 




1 02 


q 27 


8 34 





50 


D*+D- 


0.36 


0.33 


0.29 


0.2 


B^ 


D+D*- 


0.37 


0.45 


0.42 


0.9 


B^ 


D*+D- 
+ D+D*- 


0.73 


0.78 


0.71 


0.6 


B 


D^D*- 


0.54 


0.49 


0.45 


0.9 


B- 


D*°D- 


0.39 


0.36 


0.31 


0.2 


50 - 


D*+D*- 


28.78 


27.09 


25.51 


-0.1 


B 


D*^D*- 


31.37 


29.52 


27.69 


-0.1 


B^- 


D*+D*- 


1.06 


0.99 


0.91 


0.9 


B- 


D*^D*- 


1.16 


1.08 


0.98 


0.9 



change by about 10-20%, which should be the same order as annihilation and FSI effects. 
From Tables HI HTl and ITTTl there are no obvious deviations of the theoretical predictions from 
the experimental data within the present errors. It is also interesting to note that A*"^^ = 00 
is not excluded by experiments if considering the uncertainties of decay constants and from 
factors. Thus, the large N^. limit as a mechanism of factorization is not disfavored yet. 

(2) The main uncertainties of theory come from the decay constants and form factors. 
Because the decay amplitudes are proportional to decay constants, it is clear that the theo- 
retical predictions can be changed with different values of the decay constants. For instance. 
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TABLE V: Polarization fractions in various QCD approaches, where HQSi and HQSn denote the 
results with Og and ag + power corrections, respectively. 





Mode 


CQM 




LF 




HQS 


HQSi 


HQSii 




D*-D*+ 

S 


Ro = 0.523 




= 0.512 




= 0.515 


Ro = 0.517 


Ro = 0.512 






R± = 0.069 


R± 


= 0.077 


R± 


= 0.055 


= 0.070 


R± = 0.093 


B 


D*^D*- 


Ro = 0.524 


Ro 


= 0.512 


Ro 


= 0.515 


Ro = 0.517 


Ro = 0.512 






= 0.070 


R± 


= 0.078 


Ri. 


= 0.055 


= 0.070 


R± = 0.093 


B^ 


D*+D*- 


Ro = 0.547 


Ro 


= 0.535 


Ro 


= 0.538 


Ro = 0.540 


Ro = 0.535 






R± = 0.069 


R± 


= 0.077 


R± 


= 0.055 


Rj_ = 0.070 


R± = 0.092 


B- 


D*°D*- 


Ro = 0.547 


Ro 


= 0.535 


Ro 


= 0.538 


Ro = 0.541 


Ro = 0.535 






R± = 0.069 


R± 


= 0.077 


R± 


= 0.055 


R_^ = 0.070 


R± = 0.092 



the branching ratio is BR(our result) x y^^j for D^D^ . The recent experiment 

BR(i?° ^ D*+D*~) = 18.8 ± 0.9 ± 1.7 seems to favor a lower fo* ~ 0.24 than our choice of 
0.275. However, this point has to be checked by other processes. For the form factors, the 
predictions of BRs in our approach are slightly lower than those in other two approaches 
(CQM and LF). The present experiment data can not distinguish which model is more 
preferred. More precise data are necessary. Another place to test different approaches is 
through the transverse polarization R^. From Table IVl R± is predicted to be 0.07, 0.08 and 
0.09 in the CQM, LF and HQET, respectively. The larger prediction in the HQET is due 
to as corrections. Except the model-dependent calculation of power corrections in different 
approaches, one advantage of the HQET is that it permits the calculations of perturbative 
QCD corrections systematically. 

(3) The penguin effects can not be neglected in B ^ PP decays. By using the decay 
amplitudes in Appendixes, the definitions of hadronic effects in Eqs. and (fTUj) and 
the condition of ei{pi) ■ Pi = 0, we know that the effects of penguin (P) to tree (T) level, 
denoted by P/T, for PP, VP and VV modes are proportional to {af^^^ + 2a^''^^Tl)/a'f , 
{af^^^ — 2a^^^^Tl')/a'f and af^'^^ /a'f, respectively, where TZ = m'j^/[{mc + md){mb — nic)] 
and TV = m'jy/[{mc + md){mb + Trie)] and the CKM matrix elements have been canceled 
due to l^tt^i*! ~ iKfeKsl and ^ iKbKdl- The situations in the PV modes are the 
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same as those in the VV modes due to the vector meson being factorized out from the 
B —>■ P transition. Since the WCs af^'^^ and a^^'^^ have the same sign, we see clearly that 
penguin effects in the PP modes are larger than those in the VV modes; however, due to 
the cancelation between af^'^^ and a^^^\ penguin effects could be neglected in B —>■ VP 
decays. Hence, the ratios \P/T\ for PP, VP and VV{PV)) are around 15%, 0% and 4%, 
respectively. For the PP, VV{PV) modes, our predictions are consistent with the results in 
Ref 



s. Isil l32j . Note that an 4% penguin contribution was obtained for the VP modes in 
. The difference is due to that they used a lower charm quark mass (m^ = 0.95 GeV) 
than ours. For all the decay modes, the electroweak penguin contributions can be negligible 
(less than 1%). 

(4) Without FSIs, we find that the BRs in the neutral and charged modes have the 



following relationships: 

TbO Tb+ 

In addition, the decays with nonstrangeness charmed mesons are Cabibbo-suppressed com- 
pared to the decays with the Ds emission and they satisfy 



f2 

BR{B dW^W) ^ -^\^BR{B D^*^Di*^). (26) 



Clearly, if large deviations from the equalities in Eq. (|26j) are observed in experiments, they 
should arise from FSIs. Of course, if the BRs of 5° ^ L)(*)o/)(*)o ^nd 5° ^ o'f^^Di*^' 
^thO,10-V.eseen.. wn.be another MM fo.FSI.Q. 
(5) For the decay amplitude, we write 

A = T + Pe''^'^ (27) 

where T and P represent tree and penguin amplitudes, and we have chosen the convention 

such that T and P are real numbers and 6w and 6 are the CP weak and strong phases, 

respectively. From Eq. (fT^ . the CPA can be described by 

^ 2(P/r) sin ^ sin 

1 + (P/r)2 + 2(P/T) cos 6 cos Ow' ^ ^ 

According to the discussions in (1), the maximum CPAs in PP, PV and VV{VP) are 

expected to be around 26%, 0% and 8%, respectively. However, in i? ^ decays, due 
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to \ew\ = |0i|, if we take S = 90° and = 23.4", the maximum CPAs for PP and VV{VP) 
modes are 10.3% and 1.6%, respectively. Clearly, in the SM, the CPA with 0(10%) can be 
reached in i? ^ DD decays. Due to the associated CKM matrix element being Vts ~ —AX"^, 
there are no CPAs in B ^ D'^*'>Di*^ decays. In the CPA, since the strong phases mainly 
arise from the one-loop corrections which are usually small, our results on CPAs, shown in 
Table IIV[ are all at a few percent level. Therefore, if the CPAs of (9(10%) are found in 
D^D^ and B^ D^D~ decays, we can conclude the large effects of the strong phase 
are from FSIs. 

(6) As discussed before, we know that in two charmful decays the polarization fractions 
satisfy « Rq ~ R\\. The current experimental data are: Rq = 0.52 ± 0.05 P| for B^ — >■ 
D*+D*-' andi?o = 0.57±0.08±0.02, R^ = 0.19±0.08±0.01 3] and i?^ = 0.063±0.055±0.009 
jj] for B^ —>■ D*^D*~. We can see that the experimental measurements support the power- 
law relation. To estimate how large R± can be in theory, we use the relationship in Eq. ()21|) 
and the form factors in Eq. Q and we obtain 



Ri 



R 



HQS 



l + /3v + iv 

1+ Pai+ lAi 



(29) 



With the values in Eq. (j21D and R^"^^ = 0.055, we get R± ^ 10%. The detailed numerical 
values can be found in Table IVl Interestingly, for the B^ D*^D*^ decay, the estimated 

n 

result is close to the upper limit of R± = 0.063 ± 0.055 ± 0.009 observed by BABAR 4] but 
close to the lower limit of R± = 0.19±0.08±0.01 observed by BELLE 3]. We note that our 
results are different from the PQCD predictions in which R± ~ 0.06 [16|. From our results, 
we can conjecture that if large R±, say around 20%, is observed, large contributions should 
arise from FSIs. 



V. CONCLUSIONS 

We have presented a detailed study of B decaying into two charmed-mesons in the gener- 
alized factorization approach. The penguin contributions have also been taken into account. 
If the final states are both pseudoscalar mesons, the ratio of penguin and tree contributions 
is about 10% in the decay amplitude. The direct CP violating asymmetries have been esti- 
mated to be a few percent. For the -B° D*^D^, D^D*^, D*^D*^ decays, the "penguin 
pollution" is weaker than that in the D'^D~ mode. Thus, these modes provide cleaner places 
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to cross-check the value of sin2/3 measured in the 5° J/ipK decays. The weak annihila- 
tion contributions have been found to be small. We have proposed to test the annihilation 
effects in annihilation-dominated processes 

We have performed a comprehensive test of the factorization in the heavy-heavy B decays. 
The predictions of branching ratios in theory are consistent with the experimental data 
within the present level. The variations of branching ratios with the effective color number 
Nf show that the soft FSIs are not dominant. However, we cannot make the conclusion that 
they are negligible. Their effects can be of order 10-20% for branching ratios as indicated 
from the variation of N^^ . Since the soft divergences are not canceled in the non-factorizable 
corrections, this may indicate that the strong interactions at low energy either become weak 
or are suppressed by some unknown parameters (such as in the large N^. theory). If the 
factorization is still a working concept in the heavy-heavy decays, there must be some non- 
perturbative mechanisms which prefer the factorization of a large-size charmed-meson from 
an environment of "soft cloud". A relevant comment on the necessity of non-perturbative 

n 

QCD justification can be found in |84| . 

The polarization structure in the heavy-heavy decays has shown that the transverse 
perpendicular polarization fraction i?^ is the smallest while the other two are comparable 
in size. This structure follows from the QCD dynamics in the heavy quark limit. We have 
found one relation between the transverse perpendicular polarization fraction and the ratios 
of form factors, in particular V {q^) / Ai{q^) . The corrections to the heavy quark limit give 
an enhancement of R± from 0.055 to about 0.09. Since the FSIs are not significant, we do 
not expect that FSIs can change our prediction of R± substantially. If future measurements 
confirm ~ 0.2 as the recent measurement by BELLE, it will be difficult to explain within 
the HQET and the factorization hypothesis. 

In conclusion, our study has shown that the factorization works well in B meson heavy- 
heavy decays at present. More precise experimental data are desired to give a better jus- 
tification. For theory, to explain the mechanism of factorization in the heavy-heavy decays 
is of high interest. The measurement of the transverse perpendicular polarization provides 
important information on the size of the heavy quark symmetry breaking or the possibility 
of large non-factorizable effects. 
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PP DECAYS 



A{B' ^ D+D-) = v;,v;>f '''''^^ - Vt,v:, [( 



cff(c) 



X. 



{BD,Dsy 



+ 



eS{d)^{B,DDs) 



- 2a« 



efr(d)y(B,£)£),) 



)]' 



(30) 



+ {ar°' + ar^)Y, 



td 



eS{c).^{B,DsDs) 



eflf(s) ^ ^eS{c)^^{B,DsD,) 



(31) 



2^eflf(u)^(i?,DD.) 



(32) 



^Qeff(c)^(SD,D) 



2Qeff(c)^(BD,D) 



+ (af + af 



cff(c)N -^(S.DD) 



+ (af + af ^0^2 



2Qeff(d)y(B,DD) 



(33) 



+ (af")+af^yi 



cfrv(s.-DD) 



efr(u) 



+ of °' 



(34) 



A{B- ^ £>"£>-) = v;5v;>f x}""'"^ + KbKdof^^i 



cff(c) ^.{BD,D) 



'XI 



eff(c) ^{BD,Df 



- 2aT^'X. 



efr(d)^(iJ,DD) ^ eflf(d)^(B,DD) 



(35) 
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APPENDIX A: 5 ^ PV{VP) DECAYS 



+ [a 



eSv{BD,Dl) 



a 



cfr(c) ^(BD,D*) 



efr(d) ^(B,DDJ) „ efr(d)^(B,DD,') 
4 -'^l ^"6 



(Al) 



cS{d)^{B,D*Ds) 



VtbV* 



cfr(c) y{bd*,Ds) 
4 -^1 



2Qeff(d)y(s,D*i5,)^ 



(A2) 



cfr(s) 



{B,DsDl) 



(A3) 



l"3 "T "3 J-"^! 



(A4) 



* eff -^(30,01) 



cff yiB,DD*) 



+ 



eff(u)^(i3,DD*) 



a. 



efr(c) y{BD,DI) 



(A5) 



efr(c) ^{BD*,Ds) 
4 -^1 



- 2af ^*'^^)) + (af "V/^'^*^^^ - 2af '^*^^))] , (A6) 



Vcbv:, 

eflf(d) 



Vtbv:, 



a 



eff(c) y{BD,D*) 



+ («4 



efr(c) 



{B,DD*) , efr(d) 



eff(d)^(S,r'r'*) 



2ar^ ^y 



(A7) 
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D*+D- 



a. 



efr(c) ^{BD*,D) 
^1 



2 eff(c)^(i.D 



eff(d) eff(d) eff(c) 
4 "T "3 "T ^^3 



(B,D*D) 



{B,D*D) 



^^cS{d)yiB,D*D) 



(A8) 



-rp/--\ . rr/ '\ f t-> 7-v7-v*\ 



, / eff(u) I C 



(A9) 



eflf(u) 



(3,0*0) 



(AlO) 



T/ T/*^cfrx^(BD,D*) 



eS{d)^{B,DD*) 



efr(c) YiBO,D*) 



(All) 



2^eff(d)^(S,D*£)) 



(A12) 



APPENDIX "B: B^VV DECAYS 



^ D*+D*-) = K6K>f - ^tft^t 



ts 



eflf(c) y(BO\0*) 



+ 



emv{B,D*Dt) 



(Bl) 



eff(s) 



+ if' 



(B2) 
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efr(c)^(i3D*,D*) 



+ [a 



eS(n)y{B,D*D* 



(B3) 



A{B^ D*+D* 



,eff 



yiB^D-D") 



/ efr(d) cfr(d) efr(c) 

'Tyi^ ~r "3 "T "3 



I / efr(d) , . 



efr(c)^^^(B,D*D*) 



(B4) 



/^eff(u) eff(c)s^(B,D*D*) 
l"3 "T "3 J-^l 



(B5) 



eff(d)^(i?,DD) ^ eff(d)^(B,DD) 
T""4 -'1 ^"6 -'3 



(B6) 
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